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Context. The solar oxygen abundance has undergone a major downward revision in the last decade, the most noticeable one being 
^ ■ the update including 3D hydrodynamical simulations to model the solar photosphere. Up to now, such an analysis has been carried 

5— ( ' out only by one group using one radiation-hydrodynamics code. 

' Aims. We investigate the photospheric oxygen abundance considering lines from atomic transitions. We also consider the rela- 

te , tionship between the solar model used and the resulting solar oxygen abundance, to understand whether the downward abundance 

revision is specifically related to 3D hydrodynamical effects. 

Methods. We perform a new determination of the solar photospheric oxygen abundance by analysing different high-resolution 
high signal-to-noise ratio atlases of the solar flux and disc-centre intensity making use of the latest generation of C05BOLD 3D 
^ solar model atmospheres. 

OO 

Results. We find 8.73 < log(A'o/A'H) + 12 < 8.79. The lower and upper values represent extreme assumptions on the role of col- 
lisional excitation and ionisation by neutral hydrogen for the NLTE level populations of neutral oxygen. The error of our analysis 
is + (0.04± 0.03) dex, the last being related to NLTE corrections, the first error to any other effect. 3D "granulation effects" do not 
\l I play a decisive role in lowering the oxygen abundance. 

\^ , Conclusions. Our recommended value, considering our present ignorance of the role of collisions with hydrogen atoms on the 

• NLTE level populations of oxygen, is log(A'o/A?H) = 8.76 ± 0.07. The reasons which have led to lower O abundances in the past 

QQ ' are identified as (1) the lower equivalent widths adopted, and (2) the choice of neglecting collisions with hydrogen atoms in the 

^ statistical equilibrium calculations for oxygen. 

Key words. Sun: abundances - Sun: photosphere - Line: formation - hydrodynamics - convection - radiative transfer 

X 

1. Introduction In spite of its high cosmic abundance there are not many 
■ " " atomic oxygen lines which are suitable for spectroscopic abun- 
Oxygen is the most abundant element after hydrogen and he- dance determinations. Oxygen in the solar photosphere is 
Hum, and is the main product of Type II supernovae. In the mainly present in the neutral state. The permitted resonance 
studies of chemical abundances in stars, and in interstellar and and low excitation energy Unes lie in the vacuum ultraviolet, 
intergalactic matter, it is common to use the solar or solar sys- many at wavelength short-ward of the Lyman break. High- 
tem abundances as a reference. Unlike many other elements resolution spectra of the Sun at these wavelengths are not avail- 
oxygen is volatile and is therefore incompletely condensed in able. However, even if they were available, the Une blending 
meteorites, where it is present in the form of oxides. As a con- and heavy saturation of most of these Unes would make them 
sequence, the only available reference is the soUi photospheric unsuitable for abundance measurements. In the optical or near 
oxygen abundance. infrared range there are only a few O i lines which can be used, 

however, all the permitted lines are of high excitation and many 

of these are blended. There are three forbidden lines which 

Send offprint requests to: E. Caff-au ^j^^ ^^^^ ^^^^^^ ^^^^^ ^^^^ jyj^^g j^^^j^^ ^j^j^j^^ 
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however, are very weak, and they are all more or less blended 
with other atomic or molecular lines. This makes the determi- 
nation of the solar oxygen abundance a rather complex matter 
which requires special care. 

Because of its high cosmic abundance oxygen is one of the 
major contributors to the internal opacity of matter at solar-like 
chemical composition. It also plays a special role, together with 
nitrogen and carbon, in the energy production in the Sun. This 
implies that the solar oxygen abundance has a direct impact on 
the internal structure and evolutionary behaviour of the Sun. 
A solar model with reduced abun dance of CNO elements has 
lower central density and opacity jWatson & G uzik 2006) and 
thus evolves more slowly than a model with higher abundances. 
In the late evolution of the Sun, the CNO cycle becomes the 
main source of energy, and a reduced abundance of these ele- 
ments together with the lower central densities implies a higher 
central burning temperature for the Sun. 

The oscillatory properties of the Sun are strongly depen- 
dent on the metal content. In fact, it has been sugg ested to use 
helioseismic data to derive Z (lAntia & Basull2006l) . 

The solar oxygen abundance has undergone several major 
downward revi sions over the years, starting from the value of 
lLamberj (ll978h. A(0)=8.9f[ The influential review of Anders 
& Grevesse (119891) still provided a value which was near to 
this measurement (A(0)= 8.93 + 0.035), b ut nine years later it 
underwent a downward update by 0. 1 dex dGrevesse & Sauvall 
Il998i) . then was su bject to ano ther downward revision by al- 
most 0. 1 dex due to iHolweged 12001). an other downward re- 
vision bv lAllende Prieto et alj (1200 li) . and finally was lowered 
again dow n to the extremely low value of A(0)=8.66 recom- 
mended bvlAsr>lund e t al.l (l2004l). The latter result is supporte d 
by iMelendej (|2004|) and lSocas-Navarro & NortonI (l2007h . 
However, flie recent result of lAvres et al.l (l2006l A?0)=8.84). 
based on CO infrared lines, is in counter-trend. 

Th e very low abundance advocated by lAsplund et alj 
(120041) . together with the downward revision of the solar abun- 
dances of all the othe r elements by Asplund and collaborators 
dAsplund et al.ll2005l and references therein), has created seri- 
ous problems for solar models to explain the heUoseismic data 
dBahcall etaL 2005a. b: Guzik et al. 2006). 

So far, the only three analyses of the solar oxygen abun- 
dance which take into account the effects of photospheric 
inhom o ^eneities (granulation effec ts) are those of Holweger 
(2001), AUende Prieto etall ( l200lh . and that of Asplun d et alj 
C2004,) . The recommended oxygen abundances differ by almost 
0.1 dex. According to our analysis based on 3D CO^BOLD so- 
lar models, the granulation effects do not lead to a systematic 
lowering of the oxygen abundance. The 3D corrections we find 
are small and, for the majority of the lines, positive. 

We here present a new and independent determination of 
the solar oxygen abundance, with particular emphasis on a lat- 
est generation 3D CO^BOLCQ solar model, which, we believe, 
helps in understanding possible underlying systematic eiTors in 
the modelling of the solar photosphere. The line by line analy- 
sis is presented with full details in the next sections. However, 
the paper is organised in such a way that the reader not inter- 



ested in such details may directly skip to the discussion and 
conclusions in Section|6] 



2. Atomic data 

The at omic data on oxygen li nes is reviewed in detail in section 
2.1 of lAsplund et al.l (|2004 (hereafter A04). For the fines in 
common we adopt the same data and we refer the reader to that 
paper for further details. In addition to the atomic lines used by 
A04 we make use of t wo IR lines: the 1316nm line for which 
we adopt the log gf of lBiemont et al.l (Il991l) . and the 1130 nm 
line for which we use the NIST va lue of log gf = 0.076, the 
same value used bv lHolwegerl(l200lh . The adopted atomic data 
are summarised in Table [T] 

A comment on the collisional broadening con- 
stants log(76/A^H) given in Table [Tj in our 3D spectrum 
synthesis computations we payed particular attention to the 
Van der Waals broadening of the oxygen lines in order to ob- 
tain as realistic line profiles as possible. While this broadening 
process is important for the strength of strong lines it also 
plays a role for the shape of weak lines. From that perspective 
we deemed appropriate to provide the actual values which we 
employed. Whenever available we preferred Van der Waals 
broadening constants obtained from the theorv of Anstee & 
O'Mara (19953; Barklem & O'Mara (1997); Barklem et all 
(Il998al) and summarised in iBarklem et al. rdl998bl) . Hereafter, 
we shall collectively refer to the description as "ABO theory". 
ABO theory predicts a slightly different temperature depen- 
dence of the broadening effects than classical theory. It is 
customary to give the broadening constants for a Maxwellian 
(relative) velocity distribution at 10"^ K. We did not implement 
the ABO theory in our 3D spectrum synthesis code, but 
instead opted to ignore the slight differences in the temperature 
dependences and to choose a reference temperature of 5500 K 
at which we evaluated the broadening constants from ABO 
theory. This ensures that our classical implementation provided 
the ABO value at this temperature. The reference temperature 
was chosen as characteristic for the lower photospheric layers 
of the Sun. This is the reason why we give the broadening 
constants in a non-standard fashion with reference to 5500 K 
in Table □ 

Alternatively, we have computed Van der Waals broadening 
constants from the a pproximatio n implemented in the Kurucz 
routine WIDTH (see lRvanll998l) . In short, the WIDTH approx- 
imation is the following: Ce = (ae^)/h (45 - S)/Z, where S is 
the number of electrons in the ion under consideration and Z 
is the nuclear charge (for neutral oxygen 5=8 and Z=8). Then 
jJNh = ncl'\^'\ where u is the average relative velocity of 
the atom with respect to the perturber (Grav 2005, Eq. 1 1.28). 
The results are listed in Table[T] 



' A(X) = log(Nx/Nn) + 12 

^ http://www.astro.uu.se/~bf/cobold/index.html 



3. Codes and atmospheric models 

Our analysis is based on a time-dependent 3D hydrodynam- 
ical model of the solar photosphere, as well as on a number 
of standard ID hydrostatic models. The 3D model was com- 
puted with the (DO^B(p LD code (Freytag et al. 2002, 200^ 
IWedemeyer et al.ll2004l) . An overview of the numerics, aspects 
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Table 1. Atomic parameters of the oxygen lines included in 
our analysis. The table columns give the wavelength, the transi- 
tion, the oscillator strength, the excitation potential of the lower 
level, and two choices of the Van der Waals broadening con- 
stant considered: "ABO" refers to results from ABO theory, 
"WA" to the WIDTH approximation (see text for details). 



Wavelength 


Trans. 


log gf 




logijJNn) @ 5500 K 


in air [nm] 






[eV] 


ABO 


' cm^] 

WA 


630.0304 


3p-'D 


-9.717 


0.000 




-8.362 


636.3776 


3p-'D 


-10.185 


0.020 




-8.364 


615.81 




-0.296 


10.741 


-6.948 


-7.246 


777.1941 




0.369 


9.146 


-7.569 


-7.741 


777.4161 


^S-'P 


0.223 


9.146 


-7.569 


-7.741 


777.5390 




0.001 


9.146 


-7.569 


-7.741 


844.6759 


3S-3p 


0.014 


9.521 


-7.511 


-7.709 


926.66 


^^P-^D 


0.825 


10.741 


-7.334 


-7.489 


1130.2377 


5P-5S 


0.076 


10.741 


-7.141 


-7.565 


1316.4 


3P-3S 


-0.033 


10.99 


-7.108 


-7.555 



of the methodo logy of its applicatio n, and basic validity tests 
can be found in lLudwig et al.l ( l2008b . 

The 3D model used has a box size of 5. 6x5. 6x2. 27 Mm^ , 
a resolution of 140 x 140 x 150 grid points, and spans a range 
in optical depth of about -6.7 < logrRoss < 5.5. The multi- 
group opacities used by CO^BOLD are based on monochro- 
matic op acities stemming from th e MAR CS stellar atmosphere 
package jGustafsson et alj|2007l I2OO 3I) provided as function 
of gas pressure and temperature with high wavelength resolu- 
tion. For the calculation of the opacities solar elemental abun- 
dances according to Grevesse & Sauval (1998) are assumed, 
with the exceptio n of CNO for which v alues close to the rec- 
ommendation of lAsplund et al.l (|2005|) are adopted (specif- 
ically, A(C)=8.41, A(N)=7.80, A(0)=8.67). The monochro- 
matic opacities are binned into repre sentative opacity groups 
jNordlun d 1982; Ludwig 1992; L udwig et al. 1994; Vogler 
l2004l) . The solar CO^BOLD model applied here uses twelve 
bins; to our knowledge this is the highest wavelength resolu- 
tion hitherto applied in a 3D stellar atmosphere model. 

The equation of state takes into consideration the ionisa- 
tion of hydrogen and helium, as well as the formation of H2 
molecules according to Saha-Boltzmann statistics. It ignores 
the ionisation of the metals, which is of minor importance 
(<1 %) for the bulk thermodynamic properties. Relevant ther- 
modynamic quantities - in particular gas pressure and temper- 
ature - are tabulated as a function of gas density and internal 
energy. To ensure accuracy in higher derivatives (like the adi- 
abatic sound speed) analytical first derivatives are tabulated as 
well. 

Note that the tabulated EOS is used only for the computa- 
tion of the bulk thermodynamic quantities (like gas pressure, 
temperature, sound speed, specific heats) needed in the hydro- 
dynamical simulations. The electron pressure is not relevant in 
this context. In contrast, for the post facto spectrum synthesis 
from selected snapshots, a precise calculation of the electron 



pressure is important. A more detailed LTE EOS, including the 
formation of and the multiple ionisation of all individual 
metals, is employed for this purpose. 

We emphasise that our finally obtained oxygen abundance 
is not used to alter the 3D model atmosphere consistently. 
Besides the substantial computational efforts we believe that 
the feedback would be very small, in particular since oxygen 
is not a significant contributor of electrons. We tested the feed- 
back in the spectral synthesis: changing the oxygen abundance 
by 0. 1 dex or changing the oscillator strength by 0. 1 dex gave 
line strengths of the 630 nm forbidden oxygen line differing by 
only 0.1 %. 

For statistically representing the solar atmosphere, we se- 
lected 19 snapshots equidistantly spaced in time covering in to- 
tal I.2h of solar time. The snapshots are separated sufficiently 
in time so that their flow patterns are largely uncorrected. They 
have been further selected to have statistical properties close to 
the properties of the whole ensemble of available snapshots. 
We restricted ourselves to a sub-set of the available 3D data 
to keep the computational burden manageable. Here, in Fig.[Tl 
we only show a comparison of the average CO^BOLD model 
with the 3D model used by A04 in their re-determination of the 
solar oxygen abundance. The CO'' BOLD model was averaged 
on surfaces of equal Rosseland optical depth. The Rosseland 
scale was subsequently transformed to an optical depth scale 
at 500 nm in the continuum to allow a direct comparison to the 
model published by A04. This is not exactly identical to an av- 
eraging on T500 nm surfaces as done by A04 but we consider the 
involved inconsistencies immaterial for the present discussion. 
For further comparison and benchm ark, we added the empirical 
Holweger-Miiller solar model (Holweger 1967; Holweger & 
Muller ,1974„ hereafter HM) to the plot. In the low photospheric 
layers (-0.7 < logr < 0), the average CO^BOLD stratification 
follows the HM model rather closely. As a consequence, the 
CO^BOLD model reproduces the observed solar centr e-to-limb 
variation in the continuum remarkably well (see Lud wig et alj 
I2008h . In this region the A04 model has a steeper temperature 
gradient and becomes gradually cooler than the HM model. 
Towards lower optical depth, the CO^BOLD model becomes 
increasingly cooler than the HM model, while the A04 model 
approaches the HM stratification due to the different mean tem- 
perature gradients of the models. In any case, both hydrody- 
namical models stay systematically cooler than the HM model. 
In the higher photosphere (log t < -2) the CO^BOLD model 
exhibits larger temperature fluctuations than the A04 model. 

It is unclear which of the 3D models provides the better 
representation of the solar photosphere. Nonetheless, we think 
the differences seen in Fig.[T]provide an estimate of the system- 
atic uncertainties which are present in the current generation of 
hydrodynamical model atmospheres. 

As reference we used also several ID LTE plane parallel 
solar model atmospheres: 

1. An ATLAS model computed by Fiorella CastellO with 
ve rsion 9 of the ATLAS code and the solar abundances 
of lAsplund et aTl (|2005|) . With respect to other ATLAS 
solar models this has been computed with an Opacity 



http://wwwuseroats.inaf.it/castelli/sun/ap00t5777g44377klasp.dat 
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Fig. 1. Comparison of the temperature structures of the 
CO^BOLD model used in this work, the 3D model of A04, and 
the Holweger-Miiller model. For the 3D models the RMS fluc- 
tuations (thin lines) around the mean (thick lines) are indicated. 

Distribution Function ex phcitly compu t ed for the solar 
chemical composition of As plund et al.l (l2005h . a micro- 
turbulent velocity of 1 kms and a mixing-length param- 
eter Q'MLT=l-25. 

2. A solar MARCS modefl, also computed for a micro- 
turbulent velo city of 1 km s~', but with q'mlt=1-5 
dGustafsson et al. 2003, 2003). 

3. The Holweger-Miiller semi-empirical solar model 
(lHolwegeilll967l:lHolweger & Miilleill 19741) . 

4. An average 3D model denoted by (3D), generated by 
LinforSCEt this provides a ID atmospheric structure by 
averaging each 3D snapshot over surfaces of equal 
(Rosseland) optical depth (cf. Fig.[T]i. 

5. A ID hydrostatic model which employs the same micro- 
physics as C05B0LD, but uses the mi xing-length theory 
to treat convection (hereafter IDlhd, see lCafl^au & Ludwigl 
l2007bl for details). We fixed the mixing-length parameter to 
Q'MLT^I O and considered two micro-turbulence velocities, 
1.0 and l.Skms"'. We use this model as a reference for 
differential ID - 3D comparison. 

The spectrum synthesis on the 3D models is done with 
LinforSD. A detailed LTE equation of state is employed to pro- 
vide the electron pressure and the number density of all indi- 
vidual atoms and ions. A limitation of the current version of 
LinforSD is that it cannot handle a large number of lines and, 
furthermore, even if it could, the computational effort would 
be prohibitive. We used Linfor3D to compute 3D as well as ID 
synthetic spectra, for direct comparison to the hydro-model. 
For the ID synthesis we also used SYNTHE (Kurucz 1993b, 
2005bl) . in its Linux version (ISbordone et al.ll2004;.Sbordone 
2005h that has no limitation in the number of fines and tfiere- 
fore allows to include a large number (~ 1000 for a range of 
1 nm) of atomic and molecular lines in order to take into ac- 

http://marcs.astro.uu.se/ 
^ http://www.aip.de/mst/Linfor3D/Iinfor_3D_manual.pdf 
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Fig. 2. The equivalent width contribution function, 
dEW/dlogr, for the 777. Inm line at disc-centre, calcu- 
lated from the 3D model (solid line), the (3D) model (triple 
dot-dashed), a nd the IDlhd model (dashed). Following 
iMagainI (Il986h . the contribution function is defined such that 
integration over log gives the line equivalent width. 



count the blending of very w eak lines in great detail. In ad- 
dition, we employed WIDTH (Kurucz 1993al l2005bl : ICasteiiil 
.2005; . Sbordone 2005 ), to derive ID abundances directly from 
measured equivalent widths. 

The global effect produced by many weak lines is similar to 
that of an additional continuous opacity. In principle it is possi- 
ble to do an appropriate "re-normalisation" of the spectrum and 
compare this to a synthetic spectrum computed without inclu- 
sion of these very weak lines. We tested this procedure with ID 
synthetic spectra and found it to be acceptable, at least as a first 
approximation, introducing errors in the derived abundance of 
less than 0.01 dex. Therefore our approach has been to use this 
technique for the 3D analysis, keeping in the computation only 
a minimal number of lines (n < 10). 

Departures from LTE are non-negligible for several permit- 
ted O I fines, most noticeably for the 777 nm triplet. At present 
our 3D-NLTE line formation code by Steffen & Cayrel (see 
ICavrel et al.ll2007h is not capable of treating oxygen, and we 
defer full 3D-NLTE computations to a future paper 

However, in the case of Oi lines ID computations should 
yield a reliable estimate of departures from LTE. In fact the 
line formation range in the solar photosphere is very similar in 
3D and IDlhd models, respectively, for the oxygen triplet (see 
Fig.|2]i, and also for the other infrared oxygen lines. This is 
not necessarily valid for other elements and other stellar mod- 
els. For instance, 3D-NLTE and ID-NLTE contribution func- 
tions of the 670.8 nm ^Li line are very different for a metal- 
p oor model at r.fF^6240K. logq=4 .0, [M / HI ^-2.0 ( Cavrel et 
al. 120071). Even if the mean vertical stratification of oxygen is 
very similar in ID and 3D, horizontal inhomogeneities could 
introduce a difference between ID and 3D NLTE corrections. 
The error due to the ID approximation needs to be quantified 
by future detailed 3D-NLTE calculations. 
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The similarity between NLTE corrections computed in ID 
and 3D is confirmed by the results in Table 2 of A04 We sus- 
pect that the two corrections would have been even closer if the 
authors had used as a background model for their ID-NLTE 
computations a model with the mean temperature structure of 
their 3D model. 

For the present paper we decided to compute ID-NLTE 
corrections based on the ( 3D) model (rather than on IDj hp), 
and to use the Kiel code (ISteenbock & Ho lwegej 0*9841) to- 
ge ther with the oxygen model a tom describe d in Paunzen et 
al. (Il999h . We also compare our NLTE coiTections with those 
computed bv iHolweger (2001.) and A04. 

4. Observational data 

There are several high resolution, high signal-to-noise ratio 
(S/N) spectra both of the solar flux and of the disc-centre in- 
tensity, and we decided to make use of all of them. The choice 
of using more than one solar atlas is due to the fact that, even 
if S/N is very high for solar spectra, the abundances obtained 
from different spectra are not always in agreement within one 
o", highlighting the existence of systematic differences which 
are not fully understood. 

1. The data we will refer to as "Kurucz Flux" spectrum is 
based on fifty solar FTS scans taken by James Brault and 
Larry Testerman at Kitt Peak between 1981 and 1984, with 
a spectral resolving power AjAAof the order of 300 000 and 
a S/N of around 3 000, both v arying from ran ge to range. 
Further details can be found in lKurucd(l2005ah . 

2. The " Neckel Flux" and "Nec kel Intensity" spectra refer 
to the 'N eckel & Labsl (11984 ') absolutely calibrated FTS 
spectra obtained at Kitt Peak in the 1980ies, covering the 
range 330 nm to 12 50nm for ful l-disc and disc-centre, re- 
spectively (see also iNeckelll 1 999h . The spectral purity (A/l) 
ranges from 0.4pm at 330 nm to 2 pm at 1250 nm (825 000 
>/l/AA>625 000 ) 

3. Delbouille et al. ( |l973h provide a disc-centre intensity 
spectrum in the range 300 nm to lOOOnm observed from 
the JungfraujocfQ- At the same web site also an infrared 
disc-centre intens i ty spe ctrum, observed from Kitt Peak by 
iDelbouille et al.l (|1981), is available; this spectrum cov- 
ers the range from lOOOnm to 5400 nm. In the following 
we shall refer to both atlases, for brevity, as "Delbouille 
Intensity". 

4. Wallace et al. (1 19961) present d isc-cen t re inte nsity FTS 
spectra obtained at Kitt Peak by iPiercd (Il964l) . We refer 
to this data as "Wallace Intensity". 

5. For the infrared lines at 1130nm and 1316nm we anal- 
ysed Kitt Peak FTS spectra from the Digital Archive of 
the National Solar Observatory as reduced by one of the 
authors (TRA). 

The light emerging from the Sun is subject to a gravita- 
tional redshift of 636 ms"'; of the atlases we employ, only the 
Kurucz Flux has been corrected for this effect. In addition, one 
should take into account the effect of the Earth's rotation at 

* http://bass2000.obspm.fr/solar_spect.php 



the time the spectrum was acquired. Again only the Kurucz 
Flux has been carefully corrected for this effect. In principle the 
wavelength precision of the FTS scans is extremely high; that 
is, wavelengths measured at one end of the scan are very close 
relative to wavelengths measured at the other end of the scan. 
In the teiTestrial reference frame the scans have high accuracy, 
since they are slaved to a laser. However, the absolute accu- 
racy depends on how well one has corrected for all the other 
effects. Then there is the additional issue of how well the labo- 
ratory wavelengths of the oxygen lines themselves are known. 
Considering in addition that the convective line shifts predicted 
by our hydrodynamical simulations are not perfectly accurate 
(we estimate the accuracy to be of the order of 50 ms"'), we de- 
cided to introduce the wavelength shift between observed and 
synthetic spectrum as free parameter in the spectral hue fitting 
procedure described below. 

Considering the six oxygen lines at 630, 636, 777 and 
844 nm, we measure the following wavelength shifts, with 
respect to the laboratory wavelengths: for the Kurucz Flux 
21 ± 264ms-' , for the Neckel Flux 388 ± 258 ms-' and for the 
Delbouille Intensity 344 + 273 m s'' . Our r esults are in good 
agreement with Alle nde Prieto & Garcia L opez (1998) data 
who present a cr of 238 m s"' with an average of 286 m s"' for 
iron lines in the Neckel Flux spectrum. 

It is beyond the purpose of the present paper to investigate 
the convective shifts of oxygen lines. What we wanted to point 
out are the systematic differences in absolute wavelengths be- 
tween the various solar atlases. This makes it mandatory to 
leave the wavelength shift as a free parameter when performing 
line profile fitting. 

5. Data analysis 

To derive oxygen abundances, we measured the equivalent 
width (EW) of the oxygen lines and, when possible, fitted the 
line profile. EWs were computed by direct numerical (trape- 
zoidal) integration or by integration of a fitted Gaussian or 
Voigt profile, using the IRAF task splot. The results are sum- 
marised in Table[3] 

All line profile fitting has been p erformed through 
minimisation using a code described in ICaffau et all (120051) . 
Besides the oxygen abundance, also the level of the continuum 
and a wavelength shift are left as free parameters in the fitting 
procedure. For ID models, the macro turbulence has to be fixed 
a priori. The whole line profile or only a selected range can be 
fitted, as was done, for example, for the 926 nm line because of 
the presence of a residual of a telluric line in one of the wings. 

A change in the Van der Waals broadening constant y^, in- 
put to the line formation code, can produce lines with different 
shapes (more or less broadened) but also of different EW. For 
instance if we change the log of the 636 nm [OI] line from 
the standard input value for LinforSD -8.145 to the value of 
-8.364 provided by the WIDTH approximation (Ryan 199^), 
the difference in EW is 0.2%. Considering the oxygen triplet, 
the EW increases by 4% to 5% (depending on the oxygen abun- 
dance) if we change the log ye from -7.57 as derived from 
ABO theory to -7.74 as obtained from the WIDTH approxi- 
mation. We adopted the following strategy: for the transitions 
for which the broadening data were available from ABO we 
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computed the logyg according to the ABO theory; when not 
available, we used the WIDTH approximation which is the one 
used also inside the SYNTHE code. We point out that, contrary 
to naive expectations, the broadening effect is not completely 
negligible for weak lines, as also pointed out by Ryan (1998). 
In the 615 nm line, for instance, a change in log ye of 0. 1 6 dex, 
from -7.164 to -7.004, produces, for an oxygen abundance of 
8.66, a difference in EW of 1%. In contrast, the [01] lines are 
insensitive to the value adopted for log ye- 

Concerning the quadratic Stark broadening, even strong 
lines such as the oxygen triplet lines are not so sensitive to 
this parameter. In fact a change by a factor 2 in the y4 pro- 
duces in the strongest line of the triplet a change of ^o ggg dex in 
the oxygen abundance according to WIDTH combined with the 
ATLAS model. Even the shape of the line is very little affected 
by this change, only the wings being appreciably different in 
the synthetic spectra. 

Since the emphasis of this work is on the 3D analysis, we 
analysed with CO^BOLD+LinforSD all the observed spectra. 
For comparison we considered the ID computations for EWs 
for all the observed spectra, but we mostly limited the ID line- 
profile fitting to the Kurucz Flux, when available, or to the 
Delbouille Intensity. 

As in ICaffau & Ludwid (l2007bl) . we define the "3D cor- 
rections" as the difference A(0)3d - A(0)idlhd' because the 
3D effect is best captured when comparing the 3D and the 
IDlhd models (which share the same input physics) and using 
the same line formation code, LinforSD. Note that the result- 
ing "3D corrections" depend on the choice of both the mixing- 
length parameter, (Tmlt, and the micro-turbulence parameter, 
^, necessary for defining the ID model. The micro-turbulence 
is important only for the stronger lines, in our case the triplet 
lines. In this study we considered two values of the micro- 
turbulence, ^=1.0 and 1.5 km s"' (see Table [3]). Apart from 
the triplet lines, the difference in abundance when consider- 
ing the two ^ is at maximum 0.02 dex. For the triplet a different 
choice in ^ is translated to a change in the ID abundance by up 
to 0.04 dex. Higher micro-turbulence produces stronger lines 
in synthetic spectra, so that the oxygen abundance obtained is 
lower. In the present work we gave preference to ^=1.0 kms '. 
Concerning the IDlhd model for the mixing-length parameter 
we took ffMLT = 1 ■ Most of the lines are formed in layers which 
are not strongly affected by convection, and hence the choice of 
the mixing-length parameter. However, this is not true for the 
615 nm line, which is formed in deep photospheric layers. For 
testing this aspect, we considered three representative lines: the 
strongest line of the triplet at 777.1 nm, the stronger forbidden 
line at 630 nm, and the line at 615 nm. We changed the mixing- 
length parameter by ^^'^ with respect to our preferred value of 
1.0. With this change in auu the change in the IDlhd oxy- 
gen abundance (for fixed equivalent width in the synthetic flux 
spectra) for the forbidden and the triplet line is small, ^oqii 
for the forbidden line, and g^g for the triplet line. In contrast, 

for the 615 nm line we find a change in abundance of I^'q^I- 
This clearly indicates that for some lines the 3D- ID corrections 
dependent on the mixing-length parameter assumed in the ID 
comparison model. 



For completeness, we shall also provide differences using 
other ID models, either with LinforSD or with SYNTHE. We 
point out, however, that these latter differences depend also 
on code specific disparities such as different opacities, parti- 
tion functions, radiative transfer schemes, etc, and hence serve 
mainly to illustrate the magnitude of systematic modelling un- 
certainties. 



5. 1 . Corrections for departure from local thermodynamic 
equilibrium 

Besides the fact that lHolwegeil (1200 ll) computed ID non-LTE 
and A04 3D-non-LTE, the differences in the two NLTE compu- 
tations is the way in which excitation and ionisation by inelas- 
tic col lisions with hydrogen atoms is accounted for. iHolwegeJ 
('2001') accounts for them using th e formahsm of Steenbock & 
Holweger (1984), which is a generalisation of that of .Drawin I 
(11969 ). A04 instead totally ignore them. 

The two reasons which favour ignoring collisions with hy- 
drogen atoms are (i) the larger cross-section for collisions with 
electrons and (ii) the fact that electrons have much higher ve- 
locities, thereby making collisions with far more frequent. 
On the other hand, the ratio of hydrogen atoms to electrons is 
often as high as 10"* or even 10^. Under such circumstances it 
appears difficult to deny the role of collisions with hydrogen 
atoms. 

In our solar model, Nh/Ns is almost flat at lO'* in the outer 
layers, it is about 0.8 x lO'* at log tross = -0.5 and then drops 
dramatically to ~ 70 at log tr = 1 . One should therefore expect 
that the role of collisions with hydrogen atoms is different for 
lines having different line formation depths. 

Since the ISteenbock & Holweged d 19841) formalism is the 
only practical way to include the effects of collisions with hy- 
drogen atoms, it has become customary in NLTE codes to 
adopt it and to introduce the scaling factor Sh- A zero value 
for this scaling factor is equivalent to ignoring the process. For 
some species, like Na and Li, for which laboratory data exist, 
it is known that the classical Drawin formula (Drawin 196^ 
lead s to an overestimate of the coUisional rates (see Asplund 
et al. I2OO4I, and references therein). Unfortunately, no labora- 
tory data exists for collisions of neutral hydrogen with oxygen. 
Inspection of the literature shows that very different choices 
fo r Sh are made for different species. For example Gehren et 
al. (I2OOII) advocate Sh = 5 for Fe i in the Sun. 

Allende Prieto et al. (120041) have used the centre-to-limb 
variations of the 777 nm O i triplet to constrain the role of col- 
lisions with hydrogen atoms. Their analysis of synthetic 
NLTE profiles computed from their 3D hydrodynamical sim- 
ulation favours Sh = L 

In several papers the Kiel group has adopted Sh 1/3 
(Solanki & SteenbockllT988l:ISturenburg & Holwege3ll990L 
[1991 7 mainly based on unpublished work of Steenbock. In the 
present paper we present NLTE corrections for Sh - 0, 1/3, 1 
and discuss the implications. 
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5.2. Error estimates 

We here consider the possible errors that affect the analysis of 
the oxygen abundance in each line. We will take into account 
for each line the following errors: "log gf error", "broadening 
theory error", "blending error", "NLTE correction error", "sta- 
tistical error". 

1. The "log gf/ error": the uncertainty of log gf translates di- 
rectly into an error of the oxygen abundance. Since Storey 
& Zeippen (2000) do not report an error estimate of log g f 
for the two forbidden lines, we considered as the error in 
log gf the differe nce between the value given by Store y 
& Zeippen (2000) and the value used b ylLamberj (Il978h 
which has been widely used prior to the IStorev & ZeippenI 
([2OOO) computations. For the other oxygen lines we used 
the error of the /-value as listed in the NIST Atomic 
Spectra DatabasqJ- 

2. The "broadening theory error": in our analysis we applied, 
whenever possible, the ABO broadening theory. However, 
for some lines the use of the WIDTH approximation would 
result in a different abundance. We therefore consider the 
difference in the derived abundance using the two differ- 
ent broadening theories as an upper limit of the uncertainty 
associated with the broadening theory itself. 

3. The "blending error": some oxygen lines are blended with 
other lines, either stellar or telluric. Uncertainties in the 
atomic data of the blending line(s) as well as uncertain- 
ties in the abundances of the involved species introduce ad- 
ditional errors which propagate to the final oxygen abun- 
dance. In some cases, the presence of a neighbouring line, 
influencing the wings of the oxygen line in question, can 
produce an uncertainty in the placement of the contin- 
uum. This is the case, e.g., for the 844 nm line. These 
two kinds of uncertainties associated with blending were 
considered together The error due to uncertain continuum 
placement, because of blends in the line wings, has been 
estimated through successive trials, identifying the "high- 
est acceptable" and "lowest acceptable" continuum level. 
There is clearly some degree of subjectivity in this esti- 
mate. However, we were unable to devise a more objective 
approach. Concerning the telluric absorption an uncertainty 
is introduced when the line is subtracted from the spectrum. 

4. The "NLTE correction error": the choice on the value to 
assign to the H-atom collision scaling factor, Sh, changes 
the NLTE correction. As the related uncertainty we use the 
difference in the (logarithmic) NLTE abundance correction 
between S=0 and S=l/3. 

5. The "statistical error": is simply the standard deviation of 
the abundances derived from the different observed spectra. 

T he error relate d to the S/N, estimated from the Cayrel for- 
mula ( ICavrel|[T988l) . is completely negligible in comparison to 
the errors listed above. The differences between different solar 
atlases cannot be due to the random noise in the spectra, which 
is very low for all the data considered here, but must be related 
to some undetected systematic effects. 
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Fig. 3. The four observed spectra we used for the analysis of 
the [O i] 630 nm line. The Kurucz Flux spectrum, KF, and the 
Neckel Flux, NF, and Neckel Intensity, N I, spectra here p lot- 
ted are the normalised spectra provided by lKuruc^ (l2005ah and 
iNeckell ([1999) respectively; the Neckel Intensity, NI, is shifted 
up by 0.1 for clarity, while the Delbouille Intensity, DI, has 
been normalised and then shifted up by 0.1. The vertical line 
indicates the laboratory wavelength of the oxygen line. We re- 
call that KF is corrected for gravitational redshift while NF is 
not, explaining the large wavelength differences between the 
two flux atlases 



The uncertainties are summarised in Table H] The column 
"total error" of this table considers the global uncertainty that 
we attribute to each line; this value is the quadratic sum of all 
the uncertainties. The error related to the NLTE computations 
is not included in Table |4]because it is considered as a system- 
atic error Details of the derived NLTE corrections are instead 
listed in Table |5] 

5.3. Detailed analysis of individual transitions 

In this Section we give details of the analysis of each of the 
considered transitions. The results are summarised in Tables [3] 
andlH Initially, we considered also the 557 nm forbidden line, 
but discarded it because both the Delbouille Intensity and the 
Kurucz Flux show an asymmetry in the li ne profile, due t o 
the presence of C2 lines, as pointed out by iLambertl (Il978l) . 
Unfortunately the molecular data for these C2 lines is still not 
of sufficient quality to allow a reliable synthesis of the profile. 

5.3.1. [01] 630 nm 

This weak line (see Fig.O arises from the ground level and is 
immune to departures from LTE, so it has always been con- 
sidered a primary indicator of the photospheric oxygen abun- 
dance. However, it is blended with a Ni i line and lies next to 
two stronger linefl whose wings can influence the EW of the 



^ ht tp://physics.nist.gov/PhysRefData/ASD/index.html, iRalchenkol 
(I2005h 



Si I line: A = 629.9599 nm, log gf= -1.50,;rio=5.948 eV; Fe i line: 
: 630.1501 nm, log gf= -0.72,;('io=3.654eV 
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Fig. 4. The contribution of Nil line (dashed line) with 
A(Ni)=6.25 is compared to the combined [Oi]+Nii blend 
(solid line) for disc-centre according to the LinforSD compu- 
tation with A(0)=8.66. 

Table 2. Ni i isotopic splitting. 



Ni I isotope 


% 


Wavelength 


log gf 


% used 




68.27 


630.0335 


-2.11 


72.34 




26.10 


630.0355 


-2.11 


27.66 




68.27 


630.0335 


-2.11 


69.05 




26.10 


630.0355 


-2.11 


26.40 




3.59 


630.0375 


-2.11 


3.63 




0.91 


630.0395 


-2.11 


0.92 



oxygen line if not properly taken into account. In Fig. |4] the 
contribution of nickel is compared to the blend [Oi]+Nii ac- 
cording to a 3D synthetic spectrum. 

The blending nickel line exhibits a complex isotopic struc- 
ture (see Table |2]i. The two main components are ^^Ni i and 
^°Ni I, with an isotopic ratio of 0.38 and A of 630.0335 nm and 
630.0355 nm, respectively (see Bensby et al. 2004). Each com- 
ponent of the nickel bl end is assumed to have log gf--2.\\ 
([Johansson et al ] l2003h . The other two even isotopic compo- 
nents, ^^Ni I and ^'^Ni I, can be inserted in the line list as well, 
even if this does not change significantly the line shape, assum- 
ing that their shift with respect to the other even components 
has the same value as the shift between ^^Ni i and^'^Ni i ( about 
0.002 nm, see'Rosber g et al.lll993l) . [Rosberg et all (Il993h stud- 
ied the isotope splitting of the even components of Niii and 
showed that the splitting is equidistant for any of the even com- 
ponents of Ni II they studied. We assume the same behaviour 
for Nil. 

We adopted A(Ni)=6.25 (iGrevesse & Sauvall [T998h and 
consider it as a fixed parameter in the oxygen abundance de- 
termination. 



The observed solar spectra were renormalised in order to 
remove the silicon and iron contributions and the EW of the 
observed spectra were computed with the program splot. The 
oxygen abundances obtained from the curve-of-growth com- 
puted with LinforSD are reported in Table|3] Note that the abun- 
dance obtained from HM-i-Linfor3D is only slightly higher by 
a; 0.015 dex. On the other hand, both ID ATLAS and MARCS 
models, in combination with the WIDTH line formation code, 
imply a lower oxygen abundance by 0.03 dex. In this respect 
both MARCS and ATLAS models behave in a way which is 
similar to the LHD model. As can be deduced from Table [3] 
the difference 3D-(3D) is considerably smaller than the dif- 
ference 3D-1Dlhd, meaning that the latter difference is domi- 
nated by the different mean temperature structures rather than 
by horizontal temperature fluctuations. 

Note that we disagree on the sign of the 3D- ID correc- 
tion for this line with what was previously found by A04 
whose results imply a 3D correction A(0)3d - A(0)idmarcs 
of -0.04 dex for the flux, while we find a value of -1-0.03 dex 
(see above). 

Our 3D correction for this line is about +0.05 dex for the 
flux and -1-0.07 dex for the intensity spectra. The influence of 
the micro-turbulence is completely negligible for this line. 

From this line, fixing A(Ni)=6.25, the 3D oxygen abun- 
dance is 8.679 + 0.014 if we consider all four spectra. If we 
consider only the two flux spectra A(0)3d = 8.681 + 0.019, 
and A(0)3D = 8.676 + 0.014 from the two disc-centre spectra. 

Fitting the line profile of the solar intensity and flux spectra 
with a CO''BOLD-i-Linfor3D grid of synthetic spectra, we ob- 
tain an oxygen abundance which is in perfect agreement with 
the results obtained from the EWs (Table |3]l within 0.01 dex 
(see Fig.|5]l. 

The ID fits, even if not able to reproduce satisfactorily the 
line shape, are in agreement with the 3D results, giving further 
support to the notion that the 3D correction is small. The line fit 
using SYNTHE slightly underestimates the oxygen abundance 
(HM model A(0)=8.67, MARCS model A(0)=8.65 ATLAS 
model A(0)=8.61), in agreement with the above result based 
on matching the EW with ID models. 

Since this oxygen line is blended with the neutral nickel 
line, a possible approach is to use a multi parameter fit, with 
oxygen and nickel abundance, wavelength shift, and continuum 
level as free parameters. We attempted this, but the result was 
not satisfactory because wavelength shift and nickel abundance 
turned out to be closely correlated parameters, thus produc- 
ing a range of degenerate solutions for the best fit. However, 
if the wavelength scale can be established by other methods, 
the degeneracy is removed and both the nicke l and the oxy gen 
abundance can be found from line fitting (see lAyresll2008L and 
Sect.|6]below). 

Since an error in the parameters of the nickel line intro- 
duces an error of the oxygen abundance, we computed two 
additional 3D grids where we varied the Nil abundance by 
+0.1 dex with respect to its fiducial value of 6.25. Noting that 
the error of th e solar nickel abundance is 0.04 dex ( Grevesse 
& Sauval 199 8j), and the error of t he log gf value of this line 
is 0.06 dex (iJohansson et al.ll2003h . a variation of ±0.1 dex in 
the strength of the nickel line takes into account the possible 
uncertainties. Our evaluation of the systematics induced in the 
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Fig. 5. The 3D fit of the [OI] 630 nm Hne. The observed spec- 
tra (diamonds, their size indicating the observational noise 
level) are over-imposed on the best fit 3D synthetic line pro- 
file (solid). From top to bottom, the spectra are: Delbouille 
Intensity (DI), Neckel Intensity (NI), Neckel Flux (NF), and 
Kurucz Flux (KF). 



oxygen abundance shows that a difference of +0.1 in the nickel 
abundance produces a difference in the oxygen abundance of 
(resulting in A(0)=8.66, 8.74 respectively) in the fitting of 
the Kurucz Flux. Similar values are obtained for the Delbouille 
Intensity (a +0.1 change in nickel implies +0.06 change in oxy- 
gen abundance). 

Our adopted oxygen abundance for this line comes from 
the average value of the four observed spectra, and it is 
A(0)=8.68 + (0.014)s,at + (0.15)syst (see Table S for details). 
The systematic error is the quadratic sum of the log gf error 
and the eiTor related to blends and uncertainties in the contin- 
uum placement. The error due to the broadening theory is com- 
pletely negligible for this line, being of the order of 10 '^dex. 
The final result is A(0)=8.68 ± 0.15. 



- 5.3.2. [01] 636 nm 



This line lies on the red wing of a Ca i auto-ionisation line (see 
Fig.|6l). Several CN lines are present in the range. For C and 
N abundances we adopted 8.39 and 7.80 respectively. Their 
contribution to the EW of the [OI] 636 line, according to our 
ATLAS+SYNTHE synthesis, is of 0.05 pm in flux and 0.03 pm 
in intensity. With a CO^BOLD+Linfor3D synthetic spectrum, 
the CN molecules happen to come out stronger by 20% with 
respect to the ID ATLAS+SYNTHE result. This implies a dif- 
ference of only 0.002 dex in the derived oxygen abundance, 
which is far from being able to explain the 0. 1 dex difference 
between the oxygen abundance provided by this line (see be- 
low) and that of the 630 nm forbidden line (see Table[3]). 

At first we tried to remove the contribution of all the blend- 
ing lines by dividing the observed spectrum by a ID syn- 
thetic spectrum computed with all relevant atomic and molec- 
ular lines (including the Ca i auto-ionisation line) except oxy- 
gen. We obtained an oxygen abundance (A(0)=8.78), incon- 
sistent with what was derived from the other forbidden line 
(A(0)=8.68). This is at variance with what found by A04, ac- 
cording to their analysis, the difference in abundance is of only 
0.02 dex. This can stem from the differences in the EWs. But 
a difference of 0.06 pm in EW for this line can hardly explain 
the difference in A(0). A differencemust rely on the spectral- 
synthesis code used in the two analysis. Next, we normalised 
the spectrum following the profile of the Cai auto-ionisation 
line. From the EW measured from this normalised spectrum we 
subtracted the contribution due to CN molecules and obtained 
the oxygen abundances reported in Table [3] Our oxygen abun- 
dance for this line relies therefore only on the EW approach. 

We point out that the three spectra observed from Kitt 
Peak (Neckel and Kurucz) imply a higher oxygen abundance 
than the one observed from the Jungfraujoch (Delbouille). 
Moreover, also for the [OI] 630 nm line there is a difference 
in between the two intensities: in the Neckel Intensity spec- 
trum the line is stronger than in the Delbouille spectrum, even 
if the effect is less evident. We considered the hypothesis that 
the difference in height above sea level of the two sites could 
have induced a higher telluric contribution in the spectra taken 
at the site of lower elevation. But since there is no evidence 
of a telluric contribution from the inspection of spectra taken 
at different air mass, we discarded this hypothesis. More likely, 
the differences in the spectra might be related to different levels 
of solar activity at the time of the respective observations (cf. 
[Livingston et al.ll2007h . 

The fact that the oxygen abundance derived from the 
two forbidden lines is clearly not in agreement is disturbing. 
Obviously, deficiencies in the solar model atmospheres cannot 
be responsible. In fact, 3D and HM show a similar discrepancy 
(see Table [3]l. Presumably, the problem is related to our im- 
perfect knowledge of the blend l ine contribu tions. A possible 
solution is offered by the work of lAyresI (l2008i) . 

Our adopted oxygen abundance for this line, taking into 
account all the measurements, is 8.78 + (0.05)stat ± (0.12)syst, 
which is A(0)=8.78 ± 0.12 (see Table|4]for details). 
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Fig. 6. The Kurucz Flux spectrum around the [OI] 636.378 nm 
line with the identification of the stronger lines. 
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Fig. 7. The 1 615 nm Mult. 10 lines in the four observed spec- 
tra used in this work. Symbols are as in Fig. [3] The three verti- 
cal lines indicate the laboratory wavelength of the lines. 



5.3.3. Oi 615 nm 

These lines of Mult. 10 lie in a very complex spectral region 
(see Fig.|7]i, close to a rather strong iron line and in the range 
of many CN lines which make the continuum placement prob- 
lematic. The EW was measured using the de-blending capabil- 
ity of splot, with which we accounted for the presence of the 
neighbouring Fei line. The EW is very sensitive to the place- 
ment of the continuum. The oxygen abundance is then obtained 
from matching this EW with the 3D curve-of-growth. 

Our adopted LTE oxygen abundance for this line is 8.63 for 
the disc-centre spectra and 8.64 for the flux spectra. For both 
disc-centre and flux we estimate a global systematic error of 
0.21 dex, while the statistical error is 0.01 for flux and 0.03 
for intensity (see Table H] for details). We neglect any possi- 
ble error related to the molecular data because we think they 
are included in the error associated with the continuum level 
placement. 



Fig. 8. The average temperature profile (left ordinate) of the 
3D model (solid line) and the temperature profile of the IDlhd 
model (dashed) are depicted as a function of Rosseland opti- 
cal depth. The EW contribution function of the disc-centre 3D 
synthetic spectrum (right ordinate) of the 615 nm line (solid), 
is compared to that of the 630 nm (dashed) and 777 nm (triple 
dot-dashed) lines. The contribution functions are scaled to 
have similar maxima and are plotted against their respective 
monochromatic continuum optical depths. 



The IDlhd results overestimate the oxygen abundance by 
about -1-0.03 dex (flux), in other words the 3D correction is neg- 
ative, -0.03 dex. Only for this line we find a negative 3D cor- 
rection, and this is a consequence of the fact that only in this 
case the effect of horizontal temperature fluctuations, 3D-(3D), 
is appreciable. This, in turn, is a consequence of the high ex- 
citation potential of this line which implies that it forms in the 
deep photosphere and is highly sensitive to temperature fluctu- 
ations. 

Using the Holweger-Miiller model, the oxygen abundance 
is even 0.14 dex higher than the 3D result. For compari- 
son, the corresponding abundance difference is 0.12 dex for 
ATLAS -H WIDTH, and 0.06 dex for MARCS -i- WIDTH. 

Figure [8] shows the EW contribution function (following 
'Magainlll986') of the 3D disc-centre synthetic spectrum for the 
615.8 nm line (solid) in comparison with the 630.0 nm (dashed) 
and 777.1 nm (triple dot-dashed) lines, together with the tem- 
perature structures of the 3D and IDlhd models. The 615.8 nm 
line is formed deeply inside the photosphere. 

For the ID synthetic spectra we have the possibility of 
including all the relevant atomic and molecular lines and 
thus perform line profile fitting. We performed this exercise 
only on the Kurucz Flux: an ATLAS-i-SYNTHE grid provides 
A(0)=8.75. This result is consistent with those based on EW 
measurements as reported above. 

The abundance derived from this line is uncertain for the 
following reasons: the continuum placement is difficult; the 
atomic parameters of the blend lines in the range are not well 
known. We note that Holweger (2001) has discarded this line 
in his analysis of the solar oxygen abundance. We decided to 
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keep it among possible oxygen abundance indicators, although 
the large associated eiTor gives it a low weight. 

5.3.4. Oi 777 nm triplet 

These three strong lines (see Fig.|9]l lie in a clean region of 
the spectrum and do not suffer any blending from other atomic 
or molecular lines, nor of any telluric absorption. The strength 
of these lines places them on the saturation part of the curve- 
of-growth which, however, is not really flat and still provides a 
reasonable abundance sensitivity for these lines. The lines have 
relatively extended wings, such that the wings of neighbouring 
lines overlap. The EW of each line is afifected by this effect 
by 1% to 2.5% (the 777.4 nm line being the most affected). 
The triplet lines arise from a high excitation level (9.146eV) 
and are very sensitive to departures from LTE. From Table |3] 
it is evident that the LTE oxygen abundance obtained from 
CO^BOLD+LinforSD for the three Unes is larger for the fluxes 
than for the intensities, which is consistent with the fact that 
NLTE corrections are larger for flux than for intensity (see 
Table|5l). 

The 3D correction (which are expected to be sim- 
ilar in LTE and NLTE) is about +0.05 dex for flux 
and +0.10 dex for intensity (assuming a micro-turbulence 
^=1.0kms '). ATLAS+WIDTH implies the same abundances 
as CO^BOLD+LinforSD, while the abundances using the 
MARCS model are 0.06 dex lower. Therefore we conclude, 
like for that of the forbidden lines, that the direction of our 
3D correction is opposite to that found by A04 using the 
MARCS solar model as ID reference: while we find A(3D)- 
A(1Dmarcs)=+0.06, they find A(3D)-A(1Dmarcs)=-0.06. 

Since the lines are on the saturation part of the curve-of- 
growth, they are sensitive to the micro-turbulence ^ adopted 
in the ID analysis. The difference due to a change from 
l.Okms 'to 1.5 km s ' produces in both (3D) and IDlhd mod- 
els a difference of -0.04 dex (for 777. 1 nm). 

The fit of the line profile is problematic. The shape of 
the line is affected by NLTE effects (A04), so even the 
CO^BOLD+LinforSD LTE synthetic spectra can not reproduce 
very well the observed shape of the line which is wider than 
the synthetic spectra. For this reason our abundance analysis 
for these lines relies on EWs and not on line profile fitting. 

Our adopted LTE oxygen abundances for the three lines of 
the triplet as obtained by matching the measured EWs with the 
3D curve-of-growth are the following (see Tables [3] and |4] for 
details). From flux spectra: 
A(O)=9.03 + (O.OO)stat ± (0.12)syst, 

A(O)=8.99 + (0.01)stat±(0.11)syst, 

A(0)=8.96 + (O.OO)stat ± (O.lO)syst. 
From disc-centre intensity spectra: 
A(O)=8.93 + (0.01)stat±(0.12)syst, 

A(0)=8.88 + (0.02)„at ± (0.1 l)syst, 

A(O)=8.88 + (0.01)„at + (0.10)syst. 

5.3.5. Oi 844 nm 

This line (see Fig.fTOli is blended mainly by an iron line 
for which atomic data are not well known. Also in this 
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Fig. 10. The 1 844 nm line in the four observed spectra used in 
this work. Symbols are as for Fig. [3] The vertical line indicated 
the laboratory wavelength position of the line. 



case we tried to divide the observed spectrum by a ID 
ATLAS +SYNTHE synthetic one. The spectrum obtained, that 
should contain oxygen only, has a shape which is not in agree- 
ment with the CO^BOLD+LinforSD synthetic line profile in the 
blue side, where the Fe i blend line lies, but the agreement is 
good on the red side. 

A de-blending EW computation is not always stable, how- 
ever our "de-blended" EW provides an oxygen abundance in 
the range 8.80 to 8.84, depending on the observed spectrum 
considered. The fact that an over-plot of the observed spec- 
trum, "as is" with the CO^BOLD+LinforSD synthetic spectra 
suggests a value of 8.86 for the oxygen abundance and that a 
fit with the CO^BOLD+LinforSD grid avoiding the blue wing 
gives A(O)=8.80, indicates that our EW is nevertheless useful 
for the abundance determination (see Table[3]l. The 3D correc- 
tion is +0.03 dex for flux and +0.07 dex for disc-centre inten- 
sity (^ = l.Okms"'). 

Using a ID ATLAS +SYNTHE grid, we fitted both the ob- 
served Kurucz Flux spectrum "as is" and the "divided-out" 
spectrum containing only the oxygen line (obtained by divid- 
ing by a synthetic spectrum of the blend lines only), and obtain 
A(0)=8.79 and 8.73, respectively. Using a ID HM+SYNTHE 
grid in the same way, we obtain instead A(0)=8.82 and 8.76, 
respectively. 

Our adopted LTE oxygen abundance from this line is: 
A(0)=8.83 + (0.00).stat + (0.33)syst from the flux spectra, and 
A(O)=8.80 ± (0.02)stat ± (0.33)syst from disc-centi-e intensity 
spectra (see Table|4]for details). 

5.3.6. Oi 926 nm 

In the Kurucz Flux spectrum, as in the Neckel Flux and 
Intensity spectra, the Oi 926 nm line is contaminated by a 
telluric line. The intensity spectra of Delbouille (see Fig.fTTTl 
and Wallace have been processed to subtract the telluric con- 
tribution. The subtraction is probably not completely success- 
ful, since the profile still shows some distortion. Nevertheless, 
the spectra can be fitted reasonably well. Although the 
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Fig. 9. The O i 777 nm triplet in the four observed spectra used in this work. Symbols are as for Fig.|3] 
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Fig. 11. The observed Delbouille (Dl) and Wallace (WI) disc- 
centre intensity spectra of the 1 926 nm hne. 



Delbouille spectrum has a higher signal-to-noise (about 5000), 
the Wallace spectrum (S/N about 500) shows less residual dis- 
tortion after the subtraction of the telluric line. For this reason 
the EW was computed only from the Wallace intensity spec- 
trum, providing an oxygen abundance of A(0)=8.75 in 3D. The 
3D correction in this case is -1-0.06 dex - l.Okms"'). 

We also did line fitting restricted to a range that excludes 
the portion of the line that is affected by the telluric absorption. 
The fit of the Delbouille spectrum with a CO^BOLD+LinforSD 
grid gives A(0)=8.67. The 3D fitting of the Wallace spectrum 
results in A(0)=8.71. Fitting the Delbouille and Wallace spec- 
tra with an ATLAS+SYNTHE grid, we obtain A(0)=8.68 and 
8.66, respectively. With the Holweger-Miiller-i-SYNTHE grid, 
we have instead A(0)=8.74 and 8.73, respectivelj0. 

For this line our adopted LTE oxygen abundance is 
A(0)=8.74 + (O.Ol)stat + (0.12)syst (see TableSfor details). 
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Fig. 12. The observed disc-centre intensity spectrum we have 
used for the Oi 1130nm line. The vertical line indicates the 
laboratory wavelength position of the hne. 



5.3.7. O I 1130 nm 

In the Kurucz Flux atlas and in both Neckel spectra, the Oi 
1130nm line is blended by a telluric absorption. The Wallace 
intensity spectrum has been cleaned from telluric absorption, 
but its S/N is only about 80, inferior to the Kitt Peak FTS scans 
we considered, and insufficient for our purposes. 

Our analysis of this line relies therefore on a special reduc- 
tion of FTS scans obtained from the NSO Digital Archive (see 
Fig.[T2l). The Oi 1130nm line, and the companion 1316nm 
multiplet described below, fall in the shoulders of very strong 
telluric absorptions, and require a careful coiTection for the 
blending. One of us (TRA) examined the NSO/FTS catalogue 
of observations and identified a suitable series of scans, taken 
on the 13th and 14th of September 1983. Each scan was a 
co-add of two interferograms totalling 5 minutes exposure, 
through a 5" diameter aperture, at disc-centre. The datasets 
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Col. (1) spectrum identification: Dl: Delbouille Intensity, Nl: Neckel Intensity, WI: Wallace Intensity, Al: Ayres Intensity, NF: Neckel Flux, KF: 
Kurucz Flux. Col. (2) wavelength of the line. Col. (3) measured equivalent width. Col. (4) oxygen abundance, A(0), according to the CO' BOLD 
3D model. Cols. (5)-(10) A(0) from ID models, odd numbered cols, correspond to a micro-turbulence ^ of 1.0 km s"', even numbered cols, to 
f = 1.5 km s"'. Col. (1 1)-(14) 3D corrections, even numbered cols, for ^ = 1.0km s"', and odd numbered cols, for 1.5 km s"', respectively. 
(1) EW of the blend [Oi]-i-Nii; the contribution to the EWs due to Ni for A(Ni)=6.25 is 0.20 and 0.18pm for flux and intensity spectra, 
respectively, according to our 3D investigation. (2) EW relative to the two blending O i lines of the multiplet. (3) EW relative to the complete 
multiplet. 



are Usted as "8309NNrO.XXX" where "NN" is the day (13 or 
14), and "XXX" is a scan number. Hereafter these will be de- 
scribed as "NN.XXX' for simplicity. Scan 13.004 was taken 
in the early morning at about 09:00 local time, at a moderate 
airmass of 1.7. Scan 13.011 was acquired close to local noon 
at a low airmass of 1.14. Scan 13.015 was taken later in the 
afternoon, again at a moderate airmass of 1.5. On the follow- 
ing day, scan 14.007 was collected at about 09:00 local time at 
an airmass of 1.7, while scan 14.010 was taken about an hour 
before local noon, at a lower airmass of 1.27. The spectral res- 



olution of these data is very high, around 2x10^, and the S/N is 
very high as well, > lO'', taking into account only the Poisson 
noise and ignoring other sources of uncertainty. Owing to the 
excellent frequency stability of the FTS scans, pairs taken at 
low and high airmass can be registered very precisely to the 
telluric absorptions, and ratioed to eliminate the solar compo- 
nent. The ratio can be modelled as an exponential attenuation 
factor, and then applied back to the observed spectra as a cor- 
rection for the telluric component. By subjecting this technique 
to various combinations of the low and high airmass scans, one 
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Table 4. Oxygen LTE abundance for each line with the erTor estimates. The second and third column are the LTE abundance 
using the CO^BOLD 3D and HM model, respectively. The last column indicates whether a flux (F) or intensity (I) measurement 
was used. 
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Table 5. NLTE corrections for all the lines, for different values of Sh. The corrections from lHolwegeJ (1200 if) and A04 are also 
given for comparison. 
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" ID-NLTE correction from lHolwegeil ( 1200 ih 



3D-NLTE correction from A04 



can average over slight variations of the atmospheric attenua- 
tion factor in the morning versus the afternoon (owing, say, to 
increases in the humidity due to solar heating), and possible 
slight changes in the solar component due to 5 minute oscil- 
lations, convection, and the like. All these effects contribute 
to the deterioration of the information content of the spectrum 
and may be effectively thought of, as "noise". This source of 
noise is in fact much larger than the pure Poisson noise in the 
spectrum. The dispersion in the derived solar profiles is a rough 



measure of the uncertainty of the telluric correction procedure, 
and of intrinsic solar variations. For the Oi 1130nm line, the 
corresponding S/N is about 15. We refer to this S/N ratio as 
"equivalent S/N". 

The line is blended with a weaker iron line at A =1 130.242 
nm, log gf--23% \ from Kurucz's Cowan-code computations, 
whose contribution is 17% (for a ATLAS-hSYNTHE synthetic 
spectrum with A(O)=8.60) to the complete EW of the blend. 
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Fig. 13. The observed disc-centre intensity spectrum we have 
used for the analysis oftheOi 1316nm line. The vertical lines 
indicate the laboratory wavelengths of the three Oi lines dis- 
cussed in the text. 



From fitting with a CO^BOLD-i-Linfor3D grid we derive 
an oxygen abundance of A(0)=8.76. The corresponding value 
obtained from fitting with ATLAS+SYNTHE is A(0)=8.75. 
The EW based abundance determination gives a fully consis- 
tent result, A(0)=8.75. The 3D coiTection is +0.04 dex = 
l.Okms-i). 

Our adopted LTE oxygen abundance for this line is: 
A(0)=8.75 + (O.Ol)stat + (0.17)syst (see Tableglfor details). 

5.3.8. Oi 1316 nm 

As shown in Fig. [13] three oxygen lines can be found 
in this Multiplet, at A 1316.3889nm, 1316.4858nm, and 
1316.5131 nm. The strongest one (1316.4858nm) is blended 
with one of the other oxygen lines (1316.5131 nm). For the 
three lines we adopt for log g f values of -0.28, - 0.250, 
-0.728, respectively, according to lBiemont et al.l (Il99lh . 

The Wallace intensity solar spectrum has a S/N ratio of 
about 700, but the spectrum is poorly rectified, making the con- 
tinuum placement very uncertain. We decided not to use this 
spectrum but rather the Kitt Peak FTS scans mentioned above 
(Sect. 15.3^ . which have a higher S/N ratio (> 10^) and are 
well rectified. For the 1316nm feature the equivalent S/N for 
the strongest of the lines, having a central depth of about 5%, 
is about 25, somewhat higher than for the 1 130nm line. From 
fitting this spectrum, we obtain A(0)=8.75 with a 3D grid, for 
all the three lines. With an ATLAS+SYNTHE grid, the fit gives 
A(0)=8.77. 

For this line our LTE adopted abundance is A(0)=8.76 + 

(0.01)s,at±(0.22)sST. 



6. Discussion and conclusions 

The best estimate of the solar oxygen abundance is obtained by 
a suitable average of the abundances derived from the individ- 
ual lines. Our adopted oxygen abundances are summarised in 
Table |4] together with associated errors, while corrections for 



departures from LTE are listed in Table|5] A different choice in 
the collision efficiency with H-atoms, quantified by the scaling 
factor Sh (see Sec. 5.1), produces a diff'erent value in the solar 
oxygen abundance. A weighted mean 



"^M S(l/o 



(1) 



jBevington & RobinsonI [r992l p. 59), using the total errors 
from Table m provides an oxygen abundance of 



A(0)= 8.73: 
A(0) = 8.76: 
A(0) = 8.79: 



:0.04 
:0.04 
:0.04 



for Sh = 
for Sh =1/3 
for Sh = 1 



(2) 



Our preference is an intermediate value, S h=1/3, so that the so- 
lar oxygen abundance we advocate is A(0)=8.76 ± 0.04 ± 0.03. 
The first error includes the statistical uncertainties as well as 
any systematic uncertainty, except the one related to the NLTE 
corrections, which are considered separately to be 0.03 dex. If 
we sum the two errors, we arrive at an oxygen abundance of 
A(0)=8.76 ± 0.07. Computing A(0) from the simple average, 
instead of taking the weighted mean, we obtain only slightly 
smaller abundances. Oxygen abundances derived from inten- 
sity spectra are systematically smaller than those obtained from 
flux spectra, but the difference is less than 0.04 dex. 

If we consider the NLTE abundance in column 6 of Table|5] 
corresponding to Sh - 1/3, excluding the 615 nm line, we find 
a line to line scatter of 0. 13 dex, not incompatible with the Icr 
error bars we attribute to the individual lines. The same range 
in A(0) is obtained for Sh = 0. The 615 nm line requires a 
lower abundance and would increase the overall scatter, but, 
as we discussed in the section devoted to this line, this line is 
particularly problematic. 

The lowest abundance after the 615 nm line is coming from 
the forbidden oxygen line at 630 nm. The relatively low oxy- 
gen abundance it implies is sensitively dependent on the abun- 
dance of nickel which we assumed to be 6.25, the presently 
recommended value. Had be adopted the equivalent width of 
4.3 pm which A04 attribute to the oxygen contribution to the 
630 nm line we would have obtained A(0)=8.78 - moving the 
result closer to the other lines and becoming fully consistent 
with forbidden oxygen line at 636 nm. 

We do not show comparisons between line profile shapes 
of synthetic and observed spectra because we cannot reproduce 
the observed features due to blends, as is the case for the 636, 
844, and 926 nm lines; the triplet is affected by NLTE correc- 
tions. Since our NLTE computation is based on the (3D) model, 
we cannot reproduce the profile of the triplet lines accurately. 
The infra-red lines, after correction for telluric absorption, have 
low S/N so that visual comparison to synthetic profiles is not 
particularly illuminating. Only for the 630 nm forbidden line 
the observed line profile shape can be compared to the synthetic 
3D spectrum, and this is shown in Fig. |5]for all four observed 
profiles. 

Our results must be compared with those of iHolwegeJ 



(1200 Ih . A04, and to the one of UyresI (l2008l) . lHolweged(l200lh 
provides A(0)= 8.736 ± 0.078, thus -0.05 dex lower than our 
weighted mean for Sh = L right within one cr of his error bar. 
The mean oxygen abundance obtained from the atomic lines by 
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A04, A(0)=8.65, and our weighted mean, for Sh — 0, differ of 
+0.08 dex. 

One of us ( Ayres 2008) has recently carried out an inde- 
pendent analysis of specifically the 630 nm Oi forbidden line, 
taking a different approach than followed in the present work. 
Here, we have allowed the wavelength scale in the vicinity of 
the 630 nm feature to float, to accommodate the Ni i blend at 
its nominal abundance. With this approach we obtain an excel- 
lent match to the observed line-shape, although Ni i accounts 
for a larger fraction of the equivalent width of the full feature 
(see Fig.m than would be inferred from the original analysis of 
the forbidden oxygen line by AHende Prieto et al. ( 2001, their 
Fig. 1). Curiously, we obtaiii nearly the same oxygen abun- 
dance as lAllende Prieto et al.l (l200ll) . even though here the O i 
component occup ies a smaller fraction of th e full equivalent 
width. At the time'Allende Prieto et al.1 (12001') performed their 
analysis the oscillator strength of the Ni transition was not pre- 
cisely known, and they left the product gf e(Ni) as a free fitting 
parameter. They obtained a contribution corresponding to a low 
nickel abundance of A(Ni)=6.05, i.e. a rather small contribu- 
tion, using todays accepted value of the oscillator strength of 
the Ni transition. It is perceivable that the difficulty in the con- 
tinuum placemen t are responsible for the disparate findings. 

Ayres ('2008') took the alternative fitting approach of fixing 
the wavelength scale and allowing the Ni i abundance to float. 
Although the absolute velocity calibrations of the solar atlases 
and specially reduced spectra we have described earlier can be 
uncertain by hundreds of m s Ayres circumvented that de- 
ficiency by adjusting an observed FTS spectrum to match 3D- 
synthesised cores of several Fe i absorption lines in the imme- 
diate vicinity of the [O i] 630 nm feature (the 3D model was a 
single CO^BOLD snapshot in this case). This procedure moves 
the blend blue-ward of the result we obtained here, leading to a 
higher fraction of the equivalent width contributed by O i, and 
lowering the fraction of Ni i (by about 30% ). The infeiTed oxy- 
gen abundance is larger than obtained here, A(0)=8.8 1+0.04. 
Ayres's approach relies on an accurate prediction of the rela- 
tive core velocities of the Fei and Oi (+ Nil) features by the 
3D model, and assumes that the Ni i abundance can be treated 
as a free parameter. 

In both cases, it is important to note that the CO^BOLD 
model provides a different conversion between fractional O i 
equivalent width and in ferred oxygen abundance th an the 3D 
simulations reported by lAUende Prieto et al. (2001) and sub- 
sequent studies also using the Stein-Nordlund class of mod- 
els. The differences might be related to the different temper- 
ature gradients of the 3D models (e.g., as exhibited by the 
ID average stratifications shown in Fig.[Tli since a main fac- 
tor governing the strength of a weak photospheric line is the 
temperature gradient in the line formation zone. Another factor 
introducin^ystematic uncertainties is the spectral synthesis. 
[Melendez (2004) applied the spectral synthesis code MOOG 
(Sneden 1973) in his determination of the solar oxygen abun- 
dance. Besides mainly targeting OH molecular lines in the in- 
frared, he recalculated the oxygen abundance from the forbid- 
den oxygen line at 630 nm based on the average 3D model of 
A04 and their equivalent width as well as oscillator strength 
obtaining A(0)=8.75. According to our calculations horizontal 
temperature fluctuations have essentially no effect on the re- 



sulting equivalent width of the line in disc integrated spectra so 
that Melendez' value can be directly compared to the value of 
8.69 quoted by A04. We are left with the uncomfortable con- 
clusion that differences in the spectral synthesis can introduce 
differences as large as 0.06 dex in the derived oxygen abun- 
dance. 

The weighted mean of the abundances derived from several 
Oi lines in the solar spectrum implies a solar oxygen abun- 
dance in the range A(0)=8.73 to A(0)=8.79, depending on 
the actual influence of collisions by neutral hydrogen atoms on 
the NLTE level populations of oxygen. The lower abundance 
is obtained by a total neglect of such collisions, the higher 
abundance is obtained by assu ming the rather large cross sec - 
tions provided by the recipe of lSteenbock & Holweged (Il984l) . 
which may be considered an upper limit to the efficiency of 
these collisions. Any treatment of neutral hydrogen collisions 
intermediate between these two extremes will provide an oxy- 
gen abundance in the above range. 

Apart from this systematic uncertainty of the NLTE cor- 
rections, we regard the determination of A(0) as very robust, 
with a total error in the weighted mean of 0.07 dex. Considering 
our present ignorance on the actual role of collisions by neu- 
tral hydrogen, we believe that it is reasonable to assume that 
the true value lies somewhere in between these two extreme 
assumptions. Therefore, assuming Sh ~1/3, our recommended 
value for the solar oxygen abundance is the mean of the two, 
A(0) = 8.76 ± 0.07. 

Our recommended oxy gen abundance is a downwa rd revi- 
sion with respect to that of iGrevesse & Sauvall (11998 !). Using 
that set of solar abundances the solar metallicity is Zq = 0.017. 
If we take that set of solar abundances and change only the oxy- 
gen abundance to our recommended value, Zq drops to 0.016. 
On the other hand our oxygen abundance is an upward revision 
relative to what was recommended by A04. By taking the so- 
lar abundances of Asplund et al. ( 2005) but our oxygen abun- 
dance, we find Zq = 0.014. Clearly the latter result is influ- 
enced by the general downward revision of all the solar abun- 
dances advocated by Asplund and collaborators. Having found 
a higher solar oxygen abundance by using CO^BOLD models, 
and having been unable to confirm the downwar d revision of 
phosphoru s (Caffau et al. 2007) and sulphur ( Caf fau & Ludwigl 
l2007bt ICaffau et aL.2007 a). it seems not unlikely that also for 
other elements a CO^BOLD based analysis would imply higher 
abundances than recommended by I Asplund et al.l (l2005l) . 

We should compare our resul t with solar oxygen abun- 
dances deduced by other means. iDelahaye & Pinsonneaultl 
(120061) obtain A(O)=8.86±0.041 +0.025(CNNe)from the seis- 
mic properties; the latte r uncertainty is retalted to the abun- 
dance of C, N, and Ne. lAntia & Basul (l2006l) using the he- 
lioseismically measured sound speed profile obtained Z - 
0.0172 ± 0.002, but the method d oes not permit to ded uce 
the abundance of a single element. 'Chaplin et al.' (' 2007h us- 
ing the frequencies of low-degree acoustic osc illations of the 
Sun o btained Z in the range 0.0187 - 0.0229. iBasu & Antial 
(12008) find Z = 0.016, again from helioseismic considerations. 
As we can see all the results obtained from the properties of 
the solar interior favour a "high" solar metallicity, and there- 
fore a "high" oxygen abund ance. In the same line are t he re- 
sults obtained in the work of iTurck-Chieze et al ] (12004 : they 
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find a better agreement when they compare the seis mic model 
(Z-O .Ol?) to the solar model with the abundances of lHolwegej 
hoOli) (Z=0.015) with respect to the model with A04 abun- 
dances (Z=0.012). 

Th e downward revision of the solar Z by lAsplund et alJ 
(l2005h is generally considered to be mainly due to the use 
of 3D m odels, as opposed to ID models used in older anal- 
ysis (e.g. lAllende Prietdl2007h . The result of our analysis re- 
veals, however, that at least for oxygen this is not the case. The 
3D corrections for the O i lines are small and for all lines, ex- 
cept for the 615.8 nm line, they are positive, indicating that the 
3D model provides a higher abundance than the correspond- 
ing IDlhd model (see Table O. The average 3D correction, if 
we consider all the oxygen lines, is -1-0.049 assuming a micro- 
turbulence ^ = l.Okms 'in the IDlhd model, and -1-0.066 
when^ = 1.5 km s'. 

Closer inspection shows that the lower oxygen abundance 
of A04 is essentially a result of using, on average, lower equiv- 
alent widths and of the choice of Sh = (larger NLTE correc- 
tions), which are both acting in the downward direction. The 
differences between the Stein & Nordlund model employed by 
A04 and our CO'' BOLD model (see Fig.[T]i are not essential 
for the resulting oxygen abundance. The fact that our oxygen 
abundance with Sh = is 0.08 dex higher is entirely due to 
our measured equivalent widths which are larger. If we use 
the CO^BOLD model and the equivalent widths of A04 we 
find A(0)=8.67 when averaging the oxygen abundances for the 
atomic lines, very close to their result. Therefore, the 3D cor- 
rections quoted by A04 are in the opposite direction with re- 
spect to ours due to the different choice of ID reference model. 

A04 argued that their suggested low O abundance would 
place the Sun on the same oxygen level as is observed in 
the local interstellar medium and near-by B stars. However, 
it is not at all obvious that the Sun's chemical composition 
should reflect its immediate surroundings since it might have 
migr ated towards l arger galactic radii from its place of forma- 
tion (IWielen et al.iLl996) of higher metallicity. Moreover, the 
Sun harbours a planetary system and it is now well established 
that planet-hosting stars have, on average, higher metallicities 
than stars with no planets. Ecuvillon et al. (2006) presented 
a uniform analysis of oxygen abundances in a volume-limited 
sample of 59 stars with no planets and compared it to a sample 
of 96 planet-hosting stars. The mean oxygen abundance of the 
sample of planet-hosting stars is 8.77, while the mean oxygen 
abundance of the sample of stars with no planets is 8.67. Thus 
our recommended solar oxygen abundance, places the Sun in 
a very popular position among planet-hosting stars in the solar 
neighbourhood. 

To improve our knowledge of the solar oxygen abundance, 
one should focus in the first place on the collisions with hydro- 
gen atoms, both theoretically and experimentally. To the read- 
ers who have been following for some time the saga of the solar 
abundances it is perhaps not unsurprising to find that the mea- 
surement of high accuracy equivalent widths is still an open 
issue. One should try to better constrain the placement of the 
continuum, and the influence of the blending lines. 
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